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Abstract 

The attenuation of phonon modes of the DNA is due to the exchange of water 
molecules adsorbed by a molecule of DNA and the bulk solvent. Using Maxwell's 
mechanism of relaxation and a simple lattice model of the DNA, we show that 
the attenuation tends to a constant value for phonon frequencies larger than the 
inverse residence time of water molecules. We come to the conclusion that in 
the THz range the attenuation could be small enough to allow the propagation 
of the phonon modes. 



The influence of hydration on the dynamics of the 
DNA. 



Tera-hertz waves have attracted considerable attention owing to their increasing ap- 
plications and possible biological effects. There is a need of careful studying the 
action of the Tera-hertz waves on molecules of the DNA. The main point is whether 
they could generate elastic vibrations, or phonon modes, strong enough to damage a 
molecule of DNA. Woolard et al, [I], [2], have reported the direct identification of the 
phonon modes of DNA with the help of Tera-hertz spectroscopy. Theoretical analysis 
of the problem runs across difficulties because the system's small size and the range of 
hyper-sound frequencies thwart the use of the conventional mechanics of condensed 
media, [3], [1]. The problems require studying the dynamics of the duplex of the 
DNA, the so-called inter-strand modes. The experimental research and the numerical 
modelling indicated that they are strongly attenuated, [5], [6], [7], [8]. There are dif- 
ferent opinions as to the extent of the damping. In papers [9], [10], [11] the authors 
put forward arguments in favor the inter-strand modes not being over-damped. The 
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recent experimental results by Woolard et al, [T], [2], indicate that the modes are 
strong enough to be observed. 



The theoretical analysis of the elastic modes of the double helix has been per- 
formed, so far, within the framework of elasticity theory, without taking into account 
viscous effects. In this paper we shall take into account the viscosity due to the hydra- 
tion sheet, and deviations of the structure of DNA from the perfect helical symmetry, 
see [12] Thus, we shall compare a molecule of DNA with a composite material, fT5]. 

We rely on Maxwell's model of viscoelastic dynamics. According to the latter one 
may cast the relation between stress, F, and strain, u, in the form 

Here A is the elastic constant characteristic of a body, and t is the characteristic 
time of relaxation. The approach was employed for studying viscoelastic dynamics 
of liquids and amorphous bodies in the THz range of frequencies, It leads to 
important conclusions as to the mechanism of dissipation, and results in a theory 
that substantially differs from the conventional Navier-Stokes hydrodynamics. The 
most spectacular evidence to the effect is the Mandelstam-Brillouin light scatter- 
ing in liquids. The Navier-Stokes hydrodynamics gives the width of a line in the 
Mandelstam-Brilloiun doublet larger than the width of the doublet itself, so that it 
should not be observed. In fact, it is, [3j. According to the theory worked out by 
Leontovich and Mandelstam, [3], the attenuation of sound ceases to depend on its 
frequency uj in the hyper-sound region, whereas in the low frequency region, where 
the Navier-Stokes theory is valid, it varies as w^. To accommodate the hydration one 
shall recall that the two nearest layers form a kind of hydration water sheet, the bulk 
of solvent being outside. Chen and Prohofsky, [TSl suggested that the water location 
have specific sites with respect to the double helix of DNA, which form a kind of spine 
in the narrow groove. Thus, the hydration water can be considered as an integral part 
of the helix, essential to its stability and dynamics, [IS]. Viscous effects depend on 
the exchange of water molecules between the hydration sheet and the bulk of solvent, 
[16]. They have a characteristic time of tens of ps. In the spirit of the composite 
theory, |13| , we may claim that the interplay of the elastic properties of molecules 
of DNA and the hydration result in a renormalization of elastic constants. Owing to 
the THz-range of frequencies we may use Maxwell's approach to viscoelasticity, [3]. 
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2 The elastic dynamics of torsional and inter-strand 
modes 



Considering the dynamics of the DNA one has to take into account: 

(1) the DNA having the two strands; 

(2) the base-pairs being hnked by the hydrogen bonds; 

(3) the hehcal symmetry. 

To estimate the attenuation discussed above we employ a one-dimensional lattice 
model of the DNA worked out in paper PT] . The model accommodates these require- 
ments, see the Appendix. We use the scheme worked out by El Hasan and Calladine, 
[T^ , for the internal geometry of the double helix of the DNA. It describes the relative 
position of one base with respect to the other in a Watson-Crick base-pair and also 
the positions of the two base-pairs. This is achieved by introducing local frames for 
the bases and the base-pairs, and translation-slides along their long axes. We shall 
describe the relative positions of the bases of a base-pair by means of the vector Y 
directed along the long axis ( y~ axis of paper [12|). Y being equal to zero when the 
base-pair is at equilibrium. The relative positions of the base-pairs are described by 
the torsional angles 0„ , which give deviations from the standard equilibrium twist of 
the double helix. Thus a twist of the DNA molecule, which does not involve inter- 
strand motion or mutual displacements of the bases inside the pairs, is determined 
by the torsional angles (pn that are the angles of rotation of the base-pairs about the 
axis of the double-helix. The motions should correspond to the relative motion of 
the bases inside the base-pairs. For each base-pair we have the reference frame in 
which z-axis corresponds to the axis of the double helix, y-axis to the long axis of 
the base-pair, x-axis perpendicular to z- and y- axes (see Fig. 1 of paper ,12J). At 
equilibrium the change in position of adjacent base-pairs is determined only by the 
twist angle fl of the double helix. We shall assume = 27r/10 as for the B-form of 
DNA. 

The above model provides a qualitative description of the dynamics of the du- 
plex. It indicates three modes of motion: a torsional one that corresponds to the 
twist angle (/)„, and two modes for the vector Yn of inter-base displacements, [IT]- 
Accommodating the relaxation effects in accord with the recipe given above (see the 
relevant calculations in the Appendix) we may cast the dispersive equations in the 
form: 

(1) for the modes of the vectors Yn we have 



.2 _ , ,2 



aq-^ySin h - -UJaq, a = 0, 1 (2) 
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(2) for the twist mode of the duplex 



2 2 • 2 a? 1 

= sm y - -^q (3) 

Here r is the relaxation time, and wy , Ug, ujt a-re the characteristic frequencies for the 
motion of the Y, the stretching of hydrogen bonds, and the twist of the duplex (see 
the Appendix). We see that the spectrum of the twist mode has a typical acoustic 
character, whereas that for the Yn has a local minimum determined by the helical 
angle, il. Thus, the spectrum of our model is in qualitative agreement with the 
conclusions of |18| . Summarizing, we may state that there are: 

(1) the internal degrees of freedom which obey the dispersive law ([2|); 

(2) the external ones generated by winding the strands round the axis of the double 
helix, and subject to equations (|3]). 

In the range of large frequencies we may assume lot ^ 1, so that the approximate 
equation 

k^-~^ (4) 
c 2tc ^ ^ 

is true and the attenuation length. A, equals 

A w Atttc (5) 

Here c is the characteristic sound velocity of elastic excitations of a molecule of DNA. 
In fact, equation (|4]) is meaningful even for values ujt > 1, for we obtain it with the 
formula of the Newtons binom for exponent 1/2 and the imaginary correction of the 
third order. Following a rule of thumb, we may employ the last equation to find an 
estimate for the attenuation of elastic modes of DNA. 



3 Conclusions 



We see that the attenuation does not depend on frequency, if the latter is large enough, 
the length of attenuation being given by equations (|4]) and ([5]). Whether it could be 
physically significant depends on numerical values involved. The velocity of elastic 
modes of the DNA is estimated as oc 10^ cm/ sec. The relaxation time r is determined 
by the exchange of adsorbed water molecules and those of the bulk solvent. One 
may suppose that the relaxation time is of the same order of magnitude, or larger, 
as the residence time of adsorbed molecules, i.e 10^ ^^sec, according to paper [T6] . 
Phonons of a wave length of the order of 10 A have frequencies in the range of a few 
THz, and the condition a; r ^ 1 is verified. At the same time we may infer from the 
equations given above that the attenuation length is of the order A oc 10~ ^cm or 
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100 A, that is of the order of the persistence length. Therefore, the relaxation effects 
are strong, but the above estimate leaves a room for the elastic modes of DNA being 
non-over-damped. 



4 Appendix 



For the convenience of the reader we shall recall the model, [17], which relies on 
the earlier paper by H.Capellmann and W.Beim, [TO]. As was indicated above, we 
use the notations of the paper, [T^j. A twist of the DNA molecule, which does not 
involve inter-strand motion or mutual displacements of the bases inside the pairs, is 
determined by the torsional angles (j)n that are the angles of rotation of the base-pairs 
about the axis of the double-helix. The twist energy of the molecule is given by the 
equation 



E 



/ • k 



in which / is the moment of inertia, and k is the twist coefficients, which for the sake 
of simplicity and taking into account the qualitative picture at which we aim, are 
assumed to be the same for all the base-pairs. Inter-strand motions should correspond 
to the relative motion of the bases inside the base-pairs, therefore the kinetic energy 
due to this degree of freedom may be cast in the form 



E 



M .2 



where M is the effective mass of a couple. 

For each base-pair we have the reference frame in which z-axis corresponds to the 
axis of the double helix, y-axis to the long axis of the base-pair, x-axis perpendicular 
to z- and y- axes (see Fig. 1 of paper fl2l). At equilibrium the change in position of 
adjacent base-pairs is determined only by the twist angle fl of the double helix. We 
shall assume il = 27r/10 as for the B-form of DNA. To determine the energy due to 
the inter-strand displacements we need to take into account the constraint imposed 
by the helical structure of our system, which can be done as follows. Let us confine 
ourself only to the torsional degrees of freedom of the double lattice and assume the 
vectors Yn being parallel to x-y plane, or two-dimensional. Consider the displacements 
y„, Yn+i determined within the frames of the two consecutive base-pairs, n, n-l-1. 
Since we must compare the two vectors in the same frame, we shall rotate the vector 
Yn+i to the frame of the n-th base pair, 

Y^^t = R-\cj,)Y^+i 
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Here R~^{(j)) is the inverse of a 2 x 2 matrix of the rotation of the n-th frame to the 
(n+l)-one. The argument is just the same as in the KirchhofF theory of elastic rod. 
It is important that the angle (j) is given by + (j)n+i — (pn- Therefore, the energy due 
to the inter-strand stress reads 



E 



M ^ 2 



K 
2^ 



fn+l 



4'n ) Yn 



+ 1 



At equilibrium the (/)„ are equal to zero. We suppose that the size of DNA molecule 
is small enough that it can be visualized as a straight double helix, that is not larger 
than the persistence length. Hence the number of base-pairs, < 150, approximately. 
Combining the formulae given above we may write down the Lagrangian function of 
the system in the form 



E 

n 

E 



K 
2^ 



4>n)Y, 



n+1 



Y„ 



(6) 



in which K, k and a are the torsional elastic constants and the inter-pairs distance, 
correspondingly. In summations given above n is the number of a site corresponding 
to the n-th base-pair, and n = 1,2^ . . . , N , N being the number of pairs in the segment 
of the DNA under consideration. The terms, e/2Y'^ accommodate the energy of the 
inter-strand separation due to the slides of the bases inside the base-pairs. It should 
be noted that the dynamical variables (j)n and Yn are of the same order of magnitude, 
that is the first. 

Wc shall introduce the relaxation following Maxwell's prescription ([T]). We may 
cast ^ into the explicit dependence of the force on the displacement 



F = Xu - 



A(l 



1^ 
Of 



For very large relaxation times we get the Hooke law, F = Xu, and for very small 
relaxation times, that is very large dissipation, the external force is extinguished, 
F — 0. Thus, the relaxation amounts to changing scalar elastic constants A for the 
operator ones, [4], 

A = A(l - T) 

Let us consider the general situation which fully describes the problem at hand. We 
have a system given by a Lagrangian function C = C(ui,U2, • . • , Un] Ui,U2, . . . , «„; ) 



6 



In our case it is the function determined by the equations for the kinetic and the 
potential energy given above. The equations of motion have the Lagrangian form 



d_C_ _ dC_ 

dt ilk duk 



In accord with Maxwell's prescription 



we set 



dFk 
dt 



dC_ 

diik 



T 




) 



- 1 



d£ 

duk 



or 



On substituting the expression for the force in the right hand sides of the Lagrangian 
equations, and obtain the equation 



from which we infer equations (|4]),([5]). 
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